Quantum Simulation of Interacting Fermion Lattice Models in Trapped Ions 
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We propose a method of simulating efficiently many-body interacting fermion lattice models 
in trapped ions, including highly nonlinear interactions in arbitrary spatial dimensions and for 
arbitrarily distant couplings. We map products of fermionic operators onto nonlocal spin operators 
and decompose the resulting dynamics in efficient steps with Trotter methods, yielding an overall 
protocol that employs only polynomial resources. The proposed scheme can be relevant in a variety 
of fields as condensed-matter or high-energy physics, where quantum simulations may solve problems 
intractable for classical computers. 
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Quantum simulations promise to revolutionize com- 
puting technologies and scientific research by allowing 
us to solve problems that are otherwise intractable with 
classical computers [TJ [2] . Several physical systems have 
been proposed in a quantum simulator: spin models [3]- 
[6] , quantum phase transitions [7 , quantum chemistry [8] , 
anyonic statistics [9], many-body systems with Rydberg 
atoms [10] , and relativistic systems [TTHT8] . Trapped ions 
offer one of the most promising platforms for quantum 
simulators [19] , due to their high controllability, efficient 
initialization and measurement [20]. In the near future, 
experiments in quantum simulations should be able to 
solve problems intractable for classical computers, turn- 
ing this technology into a remarkable tool for scientists 
and engineers. 

The numerical simulation of fermionic systems is, 
in general, a hard problem due to the huge increase 
of the Hilbert space dimension with the number of 
modes [U [21] . Using customized numerical methods as 
quantum Monte Carlo is not always possible due to the 
well-known sign problem [22j [23] . Quantum simulations 
appear as a tool that will allow us to compute the time 
evolution of free and interacting fermion lattice theo- 
ries with minimal experimental resources. This will be 
helpful in performing a wide range of condensed mat- 
ter calculations, including those related with many-body 
interactions as the Kondo [24], Fermi-Hubbard [25], or 
Frohlich [26 Hamiltonians. Furthermore, quantum sim- 
ulations will allow us to reproduce the complete dynamics 
of these systems, avoiding mean field approximations as 
Hartree-Fock to simplify nonlinear interactions [27] . 

In this Letter, we propose a method of realizing the 
quantum simulation of many-body fermionic lattice mod- 
els for N fermionic modes in trapped ions. Our method 
can be described in three steps. Firstly, we map the set of 
N fermionic modes, via the Jordan- Wigner transforma- 
tion [28], to a set of N nonlocal spin operators. The sec- 
ond step consists in decomposing the total unitary evolu- 
tion via Trotter expansion [2, 29, 30 j in terms of a product 
of exponentials associated to each nonlocal spin operator 



appearing in the Hamiltonian. Finally, we implement 
each of these exponentials in polynomial time on a set 
of N two-level trapped ions with a reduced number of 
laser pulses [31] [32] . These three steps yield an efficient 
protocol that employs only polynomial resources. Our 
method can simulate highly nonlinear and long-range in- 
teractions in arbitrary spatial dimensions, applying the 
Jordan- Wigner transformation without the usual restric- 
tion of a reduced number of neighbors, and without the 
need of auxiliary virtual Major ana fermions [33]. This 
is due to the fact that the dynamics associated with 
the nonlocal spin operators, containing a large number 
of Pauli matrices, can still be efficiently implemented. 
The proposed protocol opens the possibility to simulate 
a wide range of interesting condensed-matter and high- 
energy physics fermionic systems for a large number of 
particles. This includes the calculation of time evolutions 
and ground state computations, e.g. through adiabatic 
protocols [34]. For a number of particles above ~ 30, 
which is foreseeable in the near future, one could already 
simulate fermionic systems that are intractable for clas- 
sical computers. 

We consider the quantum simulation of the dynamics 
associated with the general Hamiltonian 
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where Ci k has to be chosen as one of the fermionic op- 
erators bi k , b\ , that obey the anticommutation rule 

{bi k , = Si k ,i k ,, N is the number of fermionic modes, 
and a is the highest order of the many-body interaction. 

Our protocol consists of three steps, gathering tech- 
niques that have not been previously considered for quan- 
tum simulators of fermionic lattice models [35] : 

i) Jordan- Wigner mapping.- This technique estab- 
lishes a correspondence between a set of fermionic op- 
erators and a set of spin operators, transforming a lo- 
cal Hamiltonian of fermions onto a nonlocal Hamilto- 
nian of spins. Only in one spatial dimension, and for 
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couplings extended to a reduced number of neighbors, 
the correspondence is from a local model to a similar 
one [28]. The operators bi h ,b\ can always be mapped 
to products of Pauli matrices, i.e., nonlocal spin op- 
erators, using the Jordan- Wigner transformation b\ = 
In In-i ® ® ® a z k _ x ... af, and b k = (^I) 1 "- 

z(j Trotter decomposition of the total Hamiltonian. - 
Our second step consists in using standard Trotter tech- 
niques to decompose the total evolution operator in terms 
of a product of evolution operators associated to each 
nonlocal spin operator. We prove below that these evo- 
lution operators can be implemented efficiently. 

Hi) Implementation of nonlocal spin operators in 
trapped ions.- The exponentials of each nonlocal spin 
operator are efficiently implement able, given that their 
exponents consist of tensor products of k Pauli matrices. 
Each of these exponentials can be implemented, for arbi- 
trary k and up to local rotations, with a M0lmer-S0rensen 
gate upon k ions, one local gate upon one of the ions, and 
the inverse M0hner-S0rensen gate [31 j. This step can be 
summarized as 

U = Z4is(-V2, 0)Z4, (0)W M s(7t/2, 0) 

= exp a* erf erf • • • a%] , (2) 

where ^ms(#?0) = exp[— i0(cos (j)S x + sin ) 2 /4], 
&x,y = ^2^=i°~i ,y an d U az ((j)) = exp erf) for odd fc, 
where 0' = for fc = 4n + 1, and (// = —0 for k = 4n — 1, 
with positive integer n. For even /c, U az ((j)) would be 
substituted by U ay (<p) — exp(i0 / cr^), where (j)' = <p for 
fc = 4n, and 0' = — <p for fc = 4n — 2, with positive inte- 
ger n. In order to obtain directly a coupling composed of 
a y matrices times a a 2 , one may apply a similar approach 
with different M0lmer-S0rensen gates according to 

W = W M s(-7r/2,7r/2)W <T (0)W M s(7r/2,7r/2) 

= exp [i(j> af erf cr| • • • cr^] , (3) 

where the local U a (<p) gate is exp erf ) for odd fc, where 
0' = (j) for & = 4n + 1, and <fi' = —(j) for k = 4n — 1, with 
positive integer n. For even fc, the local gate is exp(i(//crf ) 
where (j)' = c/> for /c = 4n— 2, and <// = — </> for /c = 4n, with 
positive integer n. Note that local rotations acting on 
each ion give rise to any tensor product of Pauli matrices 
o~^ y,z . The coupling constant in each nonlocal spin term 
of the simulated Hamiltonian will be related to (j) 
through <p = —gt, where g is a generic coupling strength 
and t is the corresponding gate time (for details, see [35 ). 

The three steps of our protocol amounts to an efficient 
method for simulating fermionic models with long-range 
couplings in arbitrary dimensions with trapped ions. 

Note that for bounded Hamiltonians, the Trotter ex- 
pansion associated with the exponential of the polyno- 
mial sum of efficiently implementable nonlocal terms is 
also efficient [2j [29j [30], i.e., it only requires polyno- 
mial resources [35]. This includes most fermionic models 



in condensed-matter and high-energy physics, some of 
which we consider below. 

Kondo model - The long debated Kondo Hamiltonian 
provides a variety of interesting features in different sys- 
tems, as the minimum in the resistivity at low tempera- 
tures [36] . With the proposed method, we can simulate 
Kondo Hamiltonians [24] . modelling the interaction of a 
Fermi sea of electrons with several magnetic impurities 
at positions Rj, 

per pp'j 

Here, a =t,i is the spin of the electron, b PyP ^(b^ is 
the annihilation(creation) operator for an electron with 
respective momentum p or p' and spin up, a^{aj) is the 
impurity raising (lowering) spin operator, e p is the en- 
ergy of the kinetic electronic Hamiltonian and J is the 
electron-impurity coupling. Notice that, e.g., the opera- 
tors bp^bp'^aj can be now mapped to a sum of products 
of Pauli matrices, leading to an efficient implementation. 

Fermi- Hubbard model- The Fermi-Hubbard Hamil- 
tonian [25] takes into account a range of effects in 
condensed-matter physics, as the Mott transition, and is 
also believed to be relevant in high-T c superconductivity. 
It takes the form 

H = wYblbi+Sa + uYb]^]^, (5) 

5i<J j 

where the first fermionic operator subindex refers to the 
lattice site and the second to the spin, w is the hopping 
energy, U is the onsite Coulomb repulsion and one usu- 
ally makes the tight-binding approximation 5 = ±1. No- 
tice that our method is general and extends the hopping 
terms to arbitrarily distant pairs of electrons. The last 
term contains the product of four fermionic operators, 
allowing for efficient implementation. 

We could as well implement the coupling of arbitrary 
products of fermionic operators [similarly to Eq. Q] to 
linear sum of bosonic operators, 

a r N -, 

n=2 L i 1 ...i n =l j 

(6) 

The bosonic operators aj can be implemented with 
the vibrational modes of the ion chain. One would 
now consider the same gate sequence as in Eq. ([2| 
but replacing U az ((j)) = exp(^ / cr 2 ) with U <7z ^ a {(j)) = 
exp[i(f) f cr z J2j 9j( a j + a ])]- The latter can be implemented 
by using a red and a blue sideband interactions for each 
of the aj modes in the context of standard trapped-ion 
technology [37] . 

Frohlich model -The Frohlich Hamiltonian [26] models 
the interaction of electrons with phonons, 
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FIG. 1. (a) Mapping of a fermionic Hamiltonian onto an 
ion string. The couplings between fermions 1 and 4 (resp., 1 
and 10) are nonlocal when applying the Jordan- Wigner trans- 
formation, (b) Efficient mapping of the tunneling coupling 
b\bio + b\ bi in trapped ions. This highly nonlocal coupling 
can be implemented with M0lmer-S0rensen gates (dark blue 
and green), local exp(i0 > 'a%) gates (red), exp[±i(7r/4) af ] 
(yellow) and exp[±i(7r/4) af ] (cyan) gates. 



in trapped ions. This highly nonlocal coupling is a 

global unitary of 2 10 x 2 10 dimensions. In the general 

case, it would require a number of elementary gates 
of 2 20 ~ 10 6 [2 j. With our method, the number of 
gates can be as small as 10 per Trotter step, consisting 
of M0lmer-S0rensen gates (dark blue and green), local 
ex.p(i(j)' o\) gates (red), exp[±i(7r/4) J2 i a\ ] (yellow) and 
exp[±i(7r/4) ^ af ] (cyan) gates (see Fig. ^p) [35]. 

Numerical simulations. - In order to compare the effi- 
ciency of the Trotter decomposition with the exact case, 
we have realized numerical simulations of the Fermi- 
Hubbard Hamiltonian, Eq. ([5|, for different levels of 
Trotter expansion and for the exact diagonalization case. 
We have considered the case of three lattice sites, with 
six modes (two spins per site), to be simulated with six 
two- level trapped ions. The resulting Hamiltonian is 
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with M(q) being the electron-phonon coupling. Here, 
bp is the electron annihilation operator that destroys an 
electron with momentum p, a q is the phonon annihilation 
operator with momentum q, ujq is the phonon frequency, 
and m is the electron mass. The last term contains the 
product bp +q b p (a q + aL g ), whose dynamics can be imple- 
mented with our technique above according to Eq. |6]). 
We can simulate this kind of Hamiltonians in order to 
recover the polaron physics, a critical open issue for the 
deep understanding of correlated electrons in solids [38 . 

One of the main appeals of our method is that the effi- 
cient encoding of fermionic models in a lattice with arbi- 
trarily long-range couplings is feasible. This also means 
that we can apply the Jordan- Wigner transformation for 
two and three spatial dimensions, not just for one, with- 
out employing additional virtual Majorana fermions [33] . 
All this is due to the fact that the fermionic operators 
are encoded in nonlocal spin operators whose dynamics 
are efficiently implementable. Thus, the mapped spin 
Hamiltonians are highly nonlocal but their evolution is 
efficiently realizable. In order to show this, we plot in Fig. 
[I] the mapping of a solid-state 3D fermionic system onto 
an ion string. As opposite to the nearest-neighbor tun- 
neling coupling between fermions 1 and 2, which is local, 
the couplings between 1 and 4, and between 1 and 10, are 
nonlocal due to the Jordan- Wigner transformation (see 
yellow lines in Fig. [IJi). Nevertheless, we can implement 
them in an efficient way. In Fig. [l]o we show the im- 
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b 2t b ll b 2i + 4t & 3t4+ & 3+). (8) 

Notice that the number of terms to be implemented scales 
linearly with the number of modes, 57V/2 — 4 (11 in this 
case, for TV = 6). At the same time, the nonlocal gates 
upon several ions are efficiently implementable with few 
lasers, such that the number of gates in each term of the 
Hamiltonian is, in the worst case, linear in the number of 
modes, and in many cases just constant. In this partic- 
ular example the total number of gates per Trotter step 
is of 33, i.e., on average 3 gates per Hamiltonian term, 
which is very efficient. This is due to the specific struc- 
ture of this Hamiltonian, that avoids the need to apply 
additional local rotations. 

In Fig. [5J we plot in (a) and (c) the average num- 
ber of excitations for mode 621 (dashed, blue) and 63^ 
(solid, red) for different parameters, showing the good 
convergence of the Trotter method. In (a) and (c), the 
lines are obtained with exact diagonalization and the dots 
with Trotter expansion. For additional figures with larger 
number of Trotter steps, see [35]. In (b) and (c) we in- 
clude the fidelity KVK^IVK^Ot)! 2 as a function of the 
number of Trotter steps tit, where \^{t)) is the state 
evolved with exact diagonalization, and ^(^t) is the 
Trotter-evolved state, for |^(0)) = |^t(0)) = fr^JO), 
and for different parameters. Notice that the fidelity of 
the Trotter expansion goes to 1 with the number of Trot- 
ter steps nx, and for tit = 10 it is 0.99 (b) or 0.97 (d). 

With current technology, more than 100 gates have 
been realized in a single experiment [39] . Indeed, without 
error correction, one would expect the realization of more 
than thousand gates in the near future [40]. This will 
allow, for tit = 10 Trotter steps, the implementation 
of hundreds of gates per step, giving us the possibility to 
simulate a wide variety of fermionic models. In a possible 
experiment, one could consider, e.g., strings of Ca + ions 



plementation of the tunneling coupling b\bio + b\ b\ = controlled with lasers. The spin degrees of freedom can 
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FIG. 2. (a) (bl±b2i)(t) (dashed, blue), and (6j t 6 3 t>W ( solid > 
red) as a function of for a number of Trotter steps nr = 
15, and (b) fidelity \{^(tF)\^(tF)T)\ 2 as a function of tit, 
for UtF — 10, where \^(t)) is the state evolved with exact 
diagonalization, and \ip(t)r} is the Trotter-evolved state, for 
|V(0)> = \M0)) = b[^ u \0), for \w\/U = 0.1. (c) (tJ^X*) 
(dashed, blue), and (63^63^) (£) (solid, red), for tit = 15, and 
(d) fidelity \(^(t F )\^(t F ) T )\ 2 , for Ut F = 2.5, where \tp(t)) is 
the state evolved with exact diagonalization, and \ip(t)T) is 
the Trotter-evolved state, for |^(0)) = |^t(0)> = &| t &|j0), 
for |iu|/C7 = 4. In (a) and (c), the lines are obtained with 
exact diagonalization and the dots with Trotter expansion. 



be encoded in long-lived electronic states of the ions [5j|6j 
12 , 14 , 19 ; 32 , 39 . Optimal state initialization via optical 
pumping and high-fidelity detection through resonance 
fluorescence can be easily performed. 

In conclusion, we have introduced a method for simu- 
lating efficiently, i.e. with polynomial resources, many- 
body fermionic lattice models in an ion string. These 
results may be relevant for quantum simulations of 
condensed-matter systems or high-energy physics in non- 
pert urbative regimes. 
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SUPPLEMENTARY MATERIAL FOR 
QUANTUM SIMULATION OF INTERACTING 
FERMION LATTICE MODELS IN TRAPPED 
IONS" 



IMPLEMENTATION OF FERMION LATTICE 
MODELS IN TRAPPED IONS: EXAMPLES 

To illustrate our method, we propose several exam- 
ples with increasing degree of complexity. In our first 
example, in order to implement the evolution associ- 
ated with Hi = g(tf j+1 bj + b^bj + i) = g(I N g) I N -i <S> 
... <8> cr+ +1 g) crJcrT g) Ij_i g) ... g) I t + I N g) I N _i g) ... g) 
cr~ +1 g)a+a|g)/j_ig)...g)ii) in terms of M0lmer-S0rensen 
gates, we first write Hi as a sum of tensor products of 
Pauli matrices. This can always be done for Hermitian 
fermionic Hamiltonians: using a ± = (a x ± ia y )/2, we 
get Hi = g(o* +1 + ia v j+1 )/2 ® <t?(o? - «x|)/2 + (<r* +1 - 
i<jj +1 )/2 ® {&* + ia])a]/2 = -(g/2)(o* +1 a* + a v j+1 ® 
(jj). Once the fermionic Hamiltonian is in the form 
of sum of products of Pauli matrices, we use Trotter 
techniques to decompose the total evolution operator in 
product of exponentials associated with these products 
of Pauli matrices. Here, each of these exponentials is 
exp(i(g/2)t<r] +1 g) crj), and exp(i(^/2)ta] / +1 <g> crj). The 
first one can be implemented using Eq. (2) in the ar- 
ticle, with Ums(— tt/2, 0), Wms(tt/2,0) acting on ions j 
and j + 1, and U a {(j)' = —<j) = —gt/2) acting on ion 
j + 1, plus a local rotation upon ion j + 1 to change 
from the z to the x basis. Equivalently, for the sec- 
ond exponential, one can proceed similarly, substituting 
the M0lmer-S0rensen gates Ums{— tt/2, 0), Wms(tt/2,0) 
by ^ms(— tt/2, 7r/2), Ums(k/2, tt/2), acting on ions j and 
j + 1, and U(j x {(j)' = (j) = gt/2) acting on ion j + 1, plus a 



local rotation upon ion j + 1 to change from the z to the 
y basis. 

A second example describes the implementation 
of the nonlinear dynamics associated with H2 = 
g(bjbjbj +1 bj+i), where j, j + 1 may represent two dif- 
ferent modes in the same site, or in different sites. The 
j index labels the fermionic mode representing any de- 
gree of freedom (e.g., spatial, momentum, or spin). In 
this case, H2 = gcr^crj g) cr^~ +1 crj +1 , which can be writ- 
ten as a sum of products of Pauli matrices, obtaining 
H 2 = g(cr z 3 + Ij) (a| +1 + J i+1 )/4 = go* a| +1 /4 + 
gcjj/4: + <7<t| +1 /4. Finally, we apply Trotter methods 
to decompose U = exp(— iH 2 t) in three exponentials: 
exp(—igt<jj g)<r| +1 /4) will be implemented, according to 
Eq. (2) in the article, with W M s(-7r/2, 0), ^ M s(tt/2,0) 
acting on ions j and j + 1, and U ay (<fi' = —(j) = gt/4) 
acting on ion j, plus a local rotation upon ion j + 1 to 
change from the x to the z basis. The other two ex- 
ponentials, exp(— igta^/A) and exp(— zg£<rJ +1 /4) can be 
implemented easily by local gates. 

As a third example, we consider the implementation of 
the evolution associated with a tunneling Hamiltonian, 

H 3 = g(b\bi + b\ bi) = g(a% g) erf g) erf <g> ... <g> cr| <g> 

g) cr| g) erf g) ... g) erg g) cr^ ), (9) 

in trapped ions (see Fig. 1 in the article). Here, modes 
1 and 10 represent nearest neighbours in a 3D lattice. 
Note that due to the standard mapping of Jordan- Wigner 
transformation in three dimensions, the effective spin 
Hamiltonian is highly nonlocal. We first decompose the 
evolution operator ex.p(—iHst) with a Trotter expansion 
in terms of the exponentials of each of the two terms 
in i^3, i.e., exp(— igtcrf (g) erf g) erf g) ... g) erf g) crf ) and 

exip(—igtai<g>a2® a 3®---® a 9® a io)- ^ n or der to obtain the 
first of these operators, we apply the M0lmer-S0rensen 
gates Ums(— tt/2, 0), ^/ms(tt/2,0) to ions 1 to 10, apply- 
ing in between of them the local U ay (4> f = —(j) = gt) gate 
upon the second ion. Notice that, despite the large num- 
ber of ions, each of these gates requires just two lasers 
globally addressed upon all ions. We consider the sec- 
ond ion in this case because, according to Eq. (2) in 
the article, the resulting nonlocal spin operator acts with 
a a z operator upon the ion which was acted upon with 
U & ((/)' = —(j) = gt). In our case, we want it to be ion 2, 
given that in exp(— igtaf g) erf g) erf g) . . . g) erf g) af ) ion 
2, and not ion 1, is acted upon by a z operator. Instead 
of ion 2, each of the ions 3-9 could also have been chosen 
for this purpose. In this way, we reduce the amount of 
local gates needed afterwards. Accordingly, we have so 
far that W M s(-7r/2,O)Z4 v (0' = -(/> = gt)U M s^M = 
ex.p(—igtaf g) erf g) erf g) ... g) erf g) crf ). To obtain the 
desired exponential, we rotate qubits 3-9 with a global 
single qubit rotation to change the x basis to the z ba- 
sis. In order to obtain the second exponential operator, a 
similar combination of gates should be employed, but in 
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this case the specific M0lmer-S0rensen gates and single 
qubit gates would be different, as seen in Ref. [2]. 



ANALYSIS OF THE TROTTER ERROR 

In this Section we show that the resources needed in 
our protocol, including number of elementary gates and 
time scaling, are polynomial on the Trotter error, the 
total time simulated, and the total size of the system, in 
terms of the number of fermionic modes. 

The Trotter expansion is a useful tool to express the 
evolution operator of a Hamiltonian that can be written 
as a sum of efficiently-implementable Hamiltonians, in 
terms of a certain product of the operators associated to 
each of these individual Hamiltonians. More specifically, 
if the Hamiltonian H can be written as a sum of m terms, 
where m is polynomial in N, H — Y^j 
standard Trotter expansion reads [3H5], 



1 Hj , then the 



-iHt 



(]J e - iH ^ nT ) nT + 0[l/n T }. (10) 



Accordingly, by making tlt very large, the error can be 
made as small as possible. 

There are more sophisticated, higher order expansions 
so called Lie- Trotter-Suzuki methods [6-8 , that have a 
better scaling of the errors. Here, we will focus on time- 
independent Hamiltonians whose evolution operator is 
expanded in terms of a /c-th order Lie- Trotter-Suzuki in- 
tegrator. We will follow the formalism and error analysis 
of Ref. [5]. In this reference it is shown that decompo- 
sitions of U = exp(— iHt), where H = YJj=i -^j? can be 
carried out in the general form 



(ii) 



i=i 



with error \\U — U\\ < e. Here, the total number of gates 
needed, N e , scales in the error e, the total evolution time 
t, and the norm of the Hamiltonian H, \\H\\, where ||.|| 
is the 2-norm, in the form [5 

N e < m5 2fc (m||i7||t) 1+1 / 2 Ve 1/2fc , (12) 

provided e < 1 < 2m5 k - 1 \\H\\t. 

Notice that in all fermionic Hamiltonians we are con- 
sidering, we have i) a polynomial number of nonlocal spin 
operators, i.e., m is polynomial in TV, the total number 
of fermionic modes, ii) each Hj is always of the form of 
a product of arbitrary number of Pauli matrices times a 



coupling hj, such that its norm \\Hj\\ = hj, given that 
the 2-norm of a product of arbitrary number of Pauli ma- 
trices is always 1. iii) the total norm of H is bounded by 
I \ H \ I < T,7=i \\Hj\\ = Ejli hj < m^ max , where hf a * is 
the maximum among all hjS. Thus, \\H\\ is polynomial 
in m, and in consecuence, also in N. 
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FIG. 3. We plot (6^624) (t) (dashed, blue), and (6j t 6 3 t>(*) 
(solid, red) as a function of Ut, for 1^(0)) = |^t(0)) = 
6i^&i^|0), \w\/U — 4, and a number of Trotter steps nr — 30 
(a), 45 (b), and 60 (c). In (a), (b), and (c) the lines are 
obtained with exact diagonalization and the dots are ob- 
tained with Trotter expansion. We also plot (d) the fidelity 
\( r ijj{tF)\ r ^{tF)T)\ 2 as a function of 71t, for IHf — 2.5, where 
\ip(t)) is the state evolved with exact diagonalization, and 
\ip(t) T ) is the Trotter-evolved state, for |^(0)) = |^r(0)) = 
61^61^ 1 0), and for \w\/U — 4. We show the numerical results 
with Trotter (dots) and a fit to the function 1 — C/rt^ (line), 
where C is a free parameter. 



Accordingly, we have shown that the scaling of the 
number of elementary gates needed in our expansion, 
is polynomial (more specifically, a power law) in e, t, 
and N , such that our method for implementing arbitrary 
fermionic Hamiltonians that occur in nature is efficient. 



FINAL REMARKS 

We plot Fig. [3^i, b,c in order to analyze the conver- 
gence of Trotter methods to the exact diagonalization 
case when increasing the number of Trotter steps tit, and 
comparing with Fig. 2 in the article (for which tit = 15). 
These three figures clearly show the fast convergence for 
a linear increase in tit- 

In Fig. we plot the fidelity \{^(t F )\^(t F ) T )\ 2 as a 
function of n^, for Utp = 2.5, where \^(t)) is the state 
evolved with exact diagonalization, and |-0(£)t) is the 
Trotter-evolved state, for |^(0)) = |^t(0)) = fr^jjO), 
and for \w\/U = 4. We show the numerical results with 
Trotter (dots) and a fit to the function 1 — C/n^ (line), 
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where C is a free parameter. This curve has a perfect 
agreement with the Trotter numerics. Thus, the error 
goes to zero polynomially in tit, as expected. 
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